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We introduce the Blume-Emery-Griffiths (BEG) model in a social networks to describe the three- 
state dynamics of opinion formation. It shows that the probability distribution function of the 
time series of opinion is a Gaussian-like distribution. We also study the response of BEG model to 
the external periodic perturbation. One can observe that both the interior thermo-noise and the 
external field result in phase transition, which is a split phenomena of the opinion distributions. It 
is opposite between the effect acted on the opinion systems of the amplitude of the external field 
and of the thermo-noise. 
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I. INTRODUCTION 



Over the last few years, the study of opinion formation in complex networks has attracted a growing amount of 
r— !■ works and becomes the major trend of sociophysics fil. Many models have been proposed, like those of Deffuant Q], 
' Galam Krause-Hegselmann (KH) and Sznajd [||. But most models in the literature consider two-state opinion 
£pH, agents, in favor (+1) or against (—1) about a certain topic. In the Galam's majority rule and the Sznajd's updating 
' rule, the interaction between the agents is randomly changed during the evolution, and the time to reach consensus 
O ! is associated with the initial traction p of +1 state. The consensus time T reaches its maximal value at p — 0.5. In 
^ ' the Sznajd model, a pair of nearest neighbors convinces its neighbors to adopt the pair opinion if and only if both 
members have the same opinion. Otherwise the pair and its neighbors do not change opinion. In the KH consensus 
model, the opinions between and 1 and a confidence bound parameter is introduced. The agent i would take the 



average opinion of all neighboring agents that are within a confidence bound during the evolution. In the Deffuant 
Jl/' model, the opinion of two randomly selected neighboring agents i and j would remain unchanged, if their opinions 
Oi and Gj differ by more than a fixed threshold parameter. Otherwise, each opinion moves into the direction of the 
i— i , other by an amount fj, x | <jj — Oj | . 

Additionally, complex networks have received much attention in recent years. Topologically, a network is consisted 
of nodes and links. The complex network models, such as the lattice network, the random network [H, 0, HI, the 
small- world network 0, [13], and the scale- free network [ll[, are studied in many branches of science. It is meaningful 
^vq ' to mention that opinion formation models are set up in complex networks. 

In the present work, we investigate the implication of a social network in a stochastic opinion formation model. We 
£NJ ' first introduce the Blume-Emery-Griffiths (BEG) model [12, EH, E3] to describe the dynamics of opinion formation, 
and the model of complex networks we used is social network which is more reality. Our simulation focuses on the 
average opinion for different situation. And we also simulated the system under the influence of external field. 

In the rest of this paper we will give a description of this dynamic model and how to generate the underlying 
networks. In Sec lIHi we show the simulation results without external filed. In Sec. HVI we present the results with the 
' influence of external field. The final section presents further discussion and conclusion. 

II. THE MODEL 

Generally speaking, social networks include some essential characteristics, such as short average path lengths^ high 



clustering, assortative mixing [15l. Ilq|. the existence of community structure, and broad degree distributions [171 1 
As a result, we use Riitta Toivonen's social network model in our present work [19(. This network is structured by two 
processes: 1) attachment to random vertices, and 2) attachment to the neighborhood of the random vertices, giving 
rise to implicit preferential attachment. These processes give rise to essential characteristics for social networks. The 
second process gives rise to assortativity, high clustering and community structure. The degree distribution is also 
determined by the number of edges generated by the second process for each random attachment. 

In this paper, the network is grown from a chain with 10 nodes. The number of initial contacts is distributed as 
p(n>mit = 1) = 0.25, p(n,i n it = 2) = 0.75, and the number of secondary contacts from each initial contact ri2nd ~ U[0, 3] 
(uniformly distributed between and 3). The total number of nodes in the social network structure is N = 10000. The 
degree distribution of simulated networks is displayed in Fig. [TJ We note that the degree distributon P(k) is a pow er 
law functional form and a peak around the degree k = 5, also that consistent with real world observations 
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FIG. 1: Degree distribution of networks with N = 10000. Result is averages over 20 simulation runs. The number of initial 
contacts is distributed as p(rii n u = 1) = 0.25, piriinU = 2) = 0.75, and the number of secondary contacts from each initial 
Contact 722nd ~ U[0, 3]. 



Now, we consider a system with N agents, which is represented by nodes on a social network. For each node, we 
consider three states which are represented by +1, 0, and —1. A practical example could be the decision to agree 
0j(t) = +1, disagree cxj(i) = —1, or neutral <Ti(t) = 0. The states are updated according to the stochastic parallel 
spin-flip dynamics defined by the transition probabilities 

p ./ /, / 1\\ cxp{~Pej [s'\a N (t)]} 

Prob cr M+ i = s \a N {t)) = ^ , r , 77777 1 

where s, s' e {+1, 0, —1}, and j3 — a/T, a represents the active degree of system, defined as a = (a%(t)). The energy 
potential [s|<Tjv(t)] is defined by 

e* HM*)] = shi Mi)) - s 2 e % (a N (t)) , (2) 
where the following local field in node i carries all information 

h N ,i(t) = jjj<7j(t), 

On At) = E*V#)- 

Here, we define coupling and are positive numbers less than or equal to 1, and with Gaussian distribution. 
h-N.i(t) represents the time dependent interaction strengths between the node i and his rii nearest neighboring nodes. 
9N,i(t) instead the strengths of feedback and T is interior thermo-noise. So the average opinion is defined by 



K*) = ^ !>(*)• (3) 



3=1 



III. SIMULATION RESULTS 



At first we investigate the time series of average opinion, as illustrated in Fig. [2ja). It shows there exists the 
fluctuation around the average opinion r = 0. In order to compare the fluctuation of different scales, the time series 
have been normalized according to 



iZ(f) = (r(f)-(r(f)> T )/*(r(t)), 
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FIG. 2: (a) Time series of average opinion with the total time steps is t = 10000, (b) the distribution functions P(R), and (c) 
the autocorrelation function c(r). The parameters used in the simulation are p(jiinit = 1) = 0.95, N = 10000, T — 1.0 and 
L = 10000. The parameter Jij and Kij are positive numbers which are not larger than 1 in whole simulations. All the results 
in this paper are obtained over 20 realizations of the social networks. 



where (r(t)) and 5(r(t)) denote the average and the standard deviation over the period considered, respectively. In 
Fig. [D^b), we present the distribution functions P(R) associated with the time series. It is clear that this function 
P(R) is a Gaussian form. 

We calculate the autocorrelation function c(r) of our model. For a time series of L samples, r(t) for t — 1, 2, . . . , L, 
c(r) is defined by 

C(T) = Sf = 7W)-^(* + r)-f) 
Ef=T(r(i)-^) 2 

where r is the time delay and f represents the average over the period under consideration. Fig. [21c) shows the result 
of autocorrelation function of our model. It is found that c(r) decreases rapidly in very small rang of r. It means the 
system has short-time memory effects. As is now well known, the stock market has nontrivial memory effects [22| . 
For example, the autocorrelation funciton of Dow Jones (DJ), also in the small rang of r, decreases rapidly from 1 to 
0. From this point, perhaps our model is helpful to understand the financial markets. 

IV. THE INFLUENCE OF EXTERNAL FIELD 

In order to explore what phenomena maybe happen to system under the influence of external field. We add a period 
external field to the energy potential e^, 

ti [s'\a N (t)} = -shi (ctjvO)) - s 2 9 l (crjv(f)) - s [Acos(ujt + tp)] , (5) 

where A is the amplitude of period external field, to is frequency and ip denotes the initial phase of external field. 
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FIG. 3: Time series of the average opinion with different values of amplitude A = 0.08, 0.12, 0.16, 0.22, 0.28, 0.32. Parameters 
are T = 1.0, ui = n/3, and ip = 0. 

We investigate the effect of amplitude A by fixing other parameters. In Fig. [3] we plot the time series of the average 
opinion r(t) under different values of A. It is obvious that the distribution functions have a remarkable change with 
increasing A. With increasing strength of external field, the average opinion comes into several discrete parts. For 
small amplitude A = 0.02, P(R) is still a Gaussian form. When A = 0.08, it begins to appear two fluctuation around 
nonzero symmetric values of average opinions. Then, four nonzero average opinions appear at A = 0.16. Note that 
the intervals among the discrete average opinions increase with increase in the strength A of external fields. Fig. [3] 
gives the process from two wave crests to four independent parts. And the average opinion of the whole system will 
jump from one part to the other parts at all times. 

In Fig. |U we present the distribution function P(R) of the average opinion. Again, it is easy to verify that the 
average opinions oscillate among serval separate symmetric nonzero values under the external periodic driving force 
[see Fig. Ufa)]. A similar oscillation behavior is observed for simulation on the influence of the frequency u which is 
shown in Fig. SJb). Noted that P(R) for the frequency oj = 7r/3 is same to the case for u> — 2ir/3, and the same 
distribution is observed between uj = tt/6 and uj = 5tt/6. But there are distinct difference for u> = and tu = tt/2. It 
indicates a possible period ir in the case of fixed other parameters. 

Fig. [5]shows the distribution functions P(R) of average opinion time series for different initial phases ip. For ip = 0, 
the average opinion vibrates among four symmetric nonzero values. When ip increases to 7r/2, clearly, the average 
opinion comes into a 3-value oscillation. Additionally, note that the distribution functions is almost same for ip = 
and ip = 7r (or ip = tt/2 and ip — 3tt/2). Again, one can conjecture P(R) is a 7r-period behavior. We also observe the 
system's average opinion time series only have two types of distribution functions in different values of initial phases 
ip. 

Another important parameter for the systems is the interior thermo-noise T. We explore its effects with (or without) 
external fields. It is found that there is not remarkable influence on the system without external field. Contrarily, in 
the case of external field, P(R) shows a similar oscillation with it in Fig.[3Ja) (see Fig. [5]). Note that their influences 
are opposite. In Fig. [SJ with increasing T the forms of P(R) transform from four-peak to two-peak gradually, and 
merge into only one-peak at last. At the same time, the average opinion r is expanded from some separate regions to 
the whole more expansive scale for larger T. 

By comparing the Fig. SJa) with the Fig. [SJ it is clear that the amplitude A and interior thermo-noise T have 
opposite effects acting on the systems. It looks like a couple of contradictory parameters, even though both lead to 
the split phenomena of the distribution of average opinion P(R) and the nonzero average R. 

It exists similar behaviors in the Lsing ferromagnetic systems. In Ising model, the order-disorder transition is a 
second order transition. It will be a non-zero magnetization ±|M sp | for a finite system. There is a nonzero probability 
for ever that the system from near +|M sp | to near — |M sp |, and vice versa [23j ] . In our model under the influence of 
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FIG. 4: (a) The distribution functions P(R) of average opinion time series under different amplitudes A. Parameters are 
T — 1.0, u> = 7r/3, and ip = 0. (b) P(R) for different frequencies ui. Parameters are A = 0.06, <p = it/2, and T = 1.0. 

external field, it is also observed the phenomena of phase transition caused by T (or by A), which is similar to the 
Ising paramagnetic-antiferromagnetic transition. 

As discussed above, the energy potential increases with increasing T, and the system's entropy becomes larger (more 
disordered). But the external field tends to restrict the disordered effects in the system and reduces the disordered 
strength into several separate regions. 

V. CONCLUSION 

In the present work we introduce Blume-Emery-Griffiths model on opinion formation with three-state. Considering 
the characters of real social systems, we construct a social network to link between agents. In this BEG model, each 
person's opinion is influenced not only by his specific local information from his neighbors but also by the average 
opinion of the whole network. 

Moreover, we focus on the behaviors of BEG systems under external perturbation. The simulation results show 
that this system is sensitive to the external field. As discussed in Sec. IIII1 the parameters in the external periodic 
perturbation, such as amplitude A, initial phase tp, and frequency u>, have obvious impacts on the opinion systems. 
Besides, the effect of the amplitude A or interior thermo-noise T is similar to the Ising paramagnetic-antiferromagnetic 
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FIG. 5: The distribution functions P(R) of average opinion time series under different initial phases (p. Parameters are A = 0.16, 
w = it/3, and T = 1.0. 




FIG. 6: The distribution functions P(R) of average opinion time series under different interior thermo-noises T. The parameters 
used in the simulation are A — 0.16, cu — n/3, and ip = 0. 
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transition, and the influence acted on systems from A and T is opposite. 
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